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The magnetic and structural properties of MnAs are studied with ab initio methods, and by 
mapping total energies onto a Heisenberg model. The stability of the different phases is found to 
depend mainly on the volume and on the amount of magnetic order, confirming previous experimen- 
tal findings and phenomenological models. It is generally found that for large lattice constants the 
ferromagnetic state is favored, whereas for small lattice constants different antiferromagnetic states 
can be stabilized. In the ferromagnetic state the structure with minimal energy is always hexagonal, 
whereas it becomes orthorhombically distorted if there is an antiferromagnetic component in the 
hexagonal plane. For the paramagnetic state the stable cell is found to be orthorhombic up to a 
critical lattice constant of about 3.7 A, above which it remains hexagonal. This leads to the second 
order structural phase transition between paramagnetic states at about 400 K, where the lattice 
parameter increases above this critical value with rising temperature due to the thermal expan- 
sion. For the paramagnetic state an analytic approximation for the magnitude of the orthorhombic 
distortion as a function of the lattice constant is given. Within the mean field approximation the 
dependence of the Curie temperature on the volume and on the orthorhombic distortion is calcu- 
lated. For orthorhombically distorted cells the Curie temperature is much smaller than for hexagonal 
cells. This is mainly due to the fact that some of the exchange coupling constants in the hexagonal 
plane become negative for distorted cells. This is also the reason for the appearance of canted spin 
structures at low temperatures and high pressures, where the cell is orthorhombic. With these re- 
sults a description of the susceptibility as function of temperature is given, where the temperature 
dependence enters via the dependence of the Curie temperature on the lattice parameters. 



I. INTRODUCTION 

MnAs is an extremely promising material for magneto- 
electronics, since it can grow epitaxially on GaAs 1 
and Si 2 forming clean and atomically sharp interfaces* 2 - 
MnAs/GaAs heterojunctions have been extensively stud- 
ied experimentally, 1,3,4 -, and spin injection from MnAs 
into GaAs has been demonstrated. 6 However one of the 
major drawbacks for the use of MnAs in devices is the fact 
that bulk MnAs has a phase transition at 318 K, where 
the magnetic state changes from ferromagnetic to para- 
magnetic. Moreover when grown on GaAs, this temper- 
ature changes depending on the growth direction. This 
is mainly attributed to the induced straini&iS The aim 
of this paper is to use ab initio density functional the- 
ory (DFT) calculations to develop a theory of the phase 
transitions of MnAs, which can be compared with exper- 
iments and with existing phenomenological models. 

First a review of the experimental properties of MnAs 
is presented, and a brief description of the existing phe- 
nomenological models is given. Then the results of our ab 
initio calculations are presented and compared to exper- 
iments and phenomenological models. Within the scope 
of a Heisenberg model the exchange coupling constants 
are calculated for different distorted unit cells, and the 
Curie temperature and its dependence on the lattice pa- 
rameters are evaluated in the mean field approximation. 
The use of the Heisenberg model also makes possible the 
prediction of the ground state volume and lattice struc- 
ture for the paramagnetic state. It will be shown that 
the phase transitions of MnAs can indeed be explained 
by ab initio calculations. In the last section our results 
are summarized and a simple model for the susceptibility 



as function of temperature is given. From these a semi- 
quantitative description of the phase diagram of MnAs 
will emerge. 

II. EXPERIMENTAL PROPERTIES AND 
EXISTING MODELS 

A. Experimental properties 

MnAs is a ferromagnetic metal at low temperature, but 
it becomes paramagnetic at T p — 318 K. At this criti- 
cal temperature the magnetic moment changes abruptly 
from a finite value to zero (figure ^ , the resistivity in- 
creases discontinuously^ the volume is reduced by 2.1%, 
and the lattice structure changes from the hexagonal 
B8i (NiAs-type) to the orthorhombic B31 (MnP-type) 
structureiiSiiiiiSiiiii 4 * The B31 structure is obtained by 
slightly distorting the B8i structure. Associated with 
this transition there is a latent heat of 7490 J /kgiil Fur- 
thermore the phase transition shows hysteresis with a 
critical temperature of 307 K upon cooling and of 318 K 
upon heating*^ All these properties clearly indicate that 
the phase transition is of first order. 

Above T p the distortion of the crystal structure re- 
duces continuously, until it vanishes at T t « 398 KiiS At 
T t MnAs undergoes another phase transition, where it 
changes back from the B31 structure to the B8i struc- 
ture. The lattice dimensions change continuously and 
there is no latent heat, but only a discontinuity in the 
heat capacity of the material. Therefore this phase tran- 
sition is of second order. For temperatures between T p 
and T t the paramagnetic susceptibility of the material 
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FIG. 1: Magnetization per Mn atom /i as a function of P (kbar) 

temperature for ferromagnetic MnAs below 318 K, and in- 
verse susceptibility xj^ 1 for paramagnetic MnAs above 318 K FIG. 2: Temperature (T) versus pres sure (P) phase diagram 
(schematically after reference of MnAs (adapted from reference |l4lp. indicating also T p and 

Tt as functions of pressure. 



has an anomalous behavior. It increases with increasing 
temperature until it reaches a maximum at T t . Above T t 
it decreases and has a Curie- Weiss behavior (figure^. 
Moreover at T t there is a lambda point in the specific 
heat 4^ Application of a magnetic field is found to trans- 
form the B31 structure back to the B8i above a critical 
field AiSiiLi2iiS Table |U gives a short summary of the de- 
scribed properties of MnAs. 

Figure|21shows the phase diagram of MnAs. If pressure 
is applied to the system, T p is lowered, while T t increases. 
Above a critical pressure of 4.6 kbar the ferromagnetic 
B8i structure becomes unstable, and the material has 
the B31 structure for all temperatures below T t . At high 
pressures and low temperatures different types of ordered 
magnetic structures are found, with a reduced saturation 
magnetic moment with respect to the zero pressure ferro- 
magnetic phase. This, together with the anomalous be- 
havior of the susceptibility between T p and T t , led Good- 
enough to the idea that the magnetic moment of the Mn 
atoms changes from a high spin state in the B8i structure 
to a low spin state in the B31 structure i&ll According to 
his idea in the B31 structure the moment of the Mn atoms 
lies between 1 and 2 [1b ■ However neutron scattering ex- 
periments on MnAs samples in the B31 structure have 
shown that the local magnetic moment of the Mn atoms 
lies between 3 and 3.4 /xb for different temperatures and 
pressures^ in clear contradiction to the assumption of a 
low spin state. Canted spin structures, similar to the he- 
limagnetic structures of MnPj2I are found at a pressure of 
4.75 kbar below 210 K, with a local magnetic moment of 
about 3 /ie- The different saturation magnetic moments 
for different pressures therefore correspond to different 
types of canted magnetic alignments. A hysteresis region 
lies between the ferromagnetic and the canted regions, 
where both the B8i and the B31 structures can be sta- 
bilized. 

When the magnetic order breaks down and the system 
becomes paramagnetic, it has the B31 structure for all 



pressures. As the temperature is further increased the 
structure of the cell continuously changes back towards 
the B8i, until at T t it has again the B8i structure, with 
8T t /dP > 0, where P is the pressure. 

The magnetocrystalline anisotropy is quite strong in 
MnAs, with the c-axis being the hard axis, so that the 
moments prefer to lie in the hexagonal planed Measure- 
ments on the magnetoelastic coupling! 2 , indicate that the 
coupling is stronger in the hexagonal plane than perpen- 
dicular to it. 

Figure |3] shows the unit cells of MnAs in the hexago- 
nal B8i and in the orthorhombic B31 crystal structures 
and it defines the unit cell vectors a, b (t>h for the B8i 
structure) and c. The B8i structure consists of stacked 
hexagonal layers of Mn and As atoms, and the unit cell 
contains two Mn and two As atoms. The B31 structure 
has twice the volume of the B8i due to symmetry lower- 
ing and contains four Mn and four As atoms. The lattice 
is nearly hexagonal, and the atoms are moved out of the 
hexagonal symmetry points along the b and c directions 
(figure [3(b) I. 

The Mn atoms are mainly displaced in the hexago- 
nal plane along the b direction, to form chains of closer 
Mn atoms, separated from the other chains by a larger 
distance (figure 0}. The As atoms are mainly displaced 
along the c axis. In each unit cell one of the planar As 
atoms is moved upwards and the other downwards with 
respect to the original position in the B8i structure, so 
as to keep the Mn-As distance nearly constant. The dis- 
placement u of the Mn atoms in the hexagonal plane lies 
between and 0.05 b (b = |b|), depending on the tem- 
perature and pressure, and the displacement v of the As 
atoms along the c-axis lies between and 0.05 c (c = |c|). 
The B8i structure is a special case of the B31 structure, 
where b — \/3 a and u = v = 0. 

Therefore we choose the unit cell vectors in such a way 
that a and c have the same direction for both the B8i 
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Crystal structure Hexagonal B81 Orthorhombic B31 Hexagonal B81 
Magnetic order Ferromagnetic Paramagnetic Paramagnetic 

Magnetic moment 3.4/ie 

Susceptibility % dx/dT > dx/dT < 

TABLE I: Some properties of MnAs at zero pressured 



and the B31 structures. In contrast the directions of the 
vectors bh ( |bh | = |a| = a) for the hexagonal cell and b 
for the orthorhombic cell are different. 

The lattice parameters of MnAs as function of temper- 
ature and magneti c field have been measured in several 
works |!0ll8l23l24| . Figure shows the change of the 
lattice parameters as a function of temperature for zero 
pressured The lattice parameters increase with temper- 
ature due to normal thermal expansion. However the in 
plane lattice parameter a decreases when the tempera- 
tures get near T p , where it jumps from 3.717 A to the 
lower value of 3.673 A. The perpendicular lattice param- 
eter c always increases continuously with temperature. 
At Tt there is an inflection in the slope of the lattice 
parameter as a function of temperature, and at about 
T s = 450 K the slope changes discontinuously. In other 
measurements the discontinuous change of the slope is 
found at about 410 KM. 

The exact temperature at which the distortion disap- 
pears is somewhat uncertain, and fluctuations may play a 
role for small distortions. The given temperature for the 
disappearing of the distortion correspond to T t = 398 
Kji2*24 however measurements for small distortions are 
difficult and such temperature can only be inferred. As 
pointed out in reference |2j| the distortion should ap- 
pear at temperatures slightly above T t . Throughout this 
work we assume that the disappearing of the distortion 
occurs at T s , which the temperature where the thermal 
expansion coefficient of MnAs changes abruptly. 



B. Review of existing models 

Until the fifties it was believed that the first order 
phase transition at T p was between a ferromagnetic and 
an antiferromagnetic state. Kittel then proposed a model 
where the distortion from the B8i to the B31 structure 
produces the change in sign of one of the exchange cou- 
pling constants, giving rise to antiferromagnetic order m& 
However there was no experimental evidence of this anti- 
ferromagnetic state and experiments demonstrated that 
the transition is instead to a paramagnetic stated 

Therefore Bean and Rodbell (BR) proposed a modi- 
fication of Kittel's theory, where the exchange interac- 
tions are ferromagnetic for both structures, but they are 
much weaker in the B31 phase^i*2i The main idea is that 



the exchange interaction decays strongly as the volume 
decreases. In order to simplify the model, the analy- 
sis was based on the extrapolated Curie temperature Tq 
only, and not on the details of the magnetic interaction 
at the atomic level. Furthermore the model neglected the 
anisotropic change of the crystal structure, and assumed 
that it is possible to describe the change from the B8i 
to the B31 structure by a change in volume only. The 
dependence of Tq on the volume V was described by the 
following equation 



where To is the Curie temperature at the volume Vo, 
which is the volume that the system would have in the 
hypothetical absence of exchange interaction. In this con- 
text this corresponds to the volume of the B31 structure 
above T p . f3 is a parameter and is determined by fitting 
the model to experiments. Within the BR model one can 
show that for certain values of (3 a first order phase tran- 
sition between a ferromagnetic and paramagnetic state 
occurs with a simultaneous change of the volume. How- 
ever the model does not explain the second order phase 
transition at T t . Further improvements such as the in- 
troduction of a term proportional to (V — Vo) 2 do not 
change the main results^ 

Later Goodenough made an attempt to explain the 
anomalous behavior of the susceptibility and the second 
order phase transition at T t by extending the BR model 
and assuming that the local magnetic moment on the Mn 
atoms depends strongly on the volumei^i^ Here the Mn 
atoms are assumed to be in a high spin state in the B8i 
structure and in a low spin state in the B31 structure^i^ 
However measurements of the local magnetic moment 
show that the change in the magnetic moment at T p is 
very small, therefore the Goodenough model is not ap- 
plicable to this phase transition. This is probably due to 
the fact that the MnAs unit cell volume in any crystalline 
structure is too big to justify a high-spin to low-spin tran- 
sition. In fact in MnAsi-^P^ the change from high-spin 
to low-spin state is observedrfl*22i2!i however the low spin 
state is found for unit cell volumes smaller than 120 A 3 . 
In contrast the unit cell of MnAs has always a volume of 
around 130 A 3 , and it is always in the high-spin state. 

In 1982 Kato et al. extended the BR model by tak- 
ing into account not only the change in volume, but also 
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FIG. 3: (a) B81 unit cell containing two Mn and two As 
atoms, (b) B31 unit cell containing four Mn and four As 
atoms. 



the change of the crystal structured Furthermore in- 
stead of just using Tq for describing the magnitude of 
the exchange interactions, they consider exchange cou- 
pling constants up to second nearest neighbors. The ob- 
tained results are similar to those of Bean and Rodbell. 
The second order phase transition at T t is explained by 
assuming that locally the structure above T t is still the 
B31, but that the distortions from the B81 are randomly 
distributed. This assumption is probably not valid, since 




FIG. 4: Two-dimensional representation of one layer of Mn 
atoms in the B31 structure, dij represent the various Mn-Mn 
distances. The first index i=2 indicates that all the Mn are 
second nearest-neighbor to each other in the B8i structure. 
The second index j =1, 2, 3 labels the three distances arising 
from the B31 distortion. 

neutron diffraction experiments clearly indicate that the 
structure is a regular B8i above T t i2i 

A Landau-type phenomenological model, where the 
free energy $ is expanded as a functio n of order pa- 
rameters, is given in references |25l32l33| and references 
therein. In reference |25| just two order parameters are 
used, the relative magnetization a (a — 1 for a ferromag- 
netic state, er = for a paramagnetic state) and the or- 
thorhombic distortion d. The following equation is used 
for the free energy: 

*(d, ct;T, H) = $ + ci (T - T (1 - <5i d 2 )) a 2 + c 2 a 4 + 
+c 3 cr 6 + c 4 (T- To) d 2 + c 5 d 4 - M aH (l - 8 2 d 2 ) , 

(2) 

where ci, C5 and Si, 62 are expansion coefficients to be 
fitted to experiment, H is an external magnetic field, Mq 
is the saturation magnetization, To is the extrapolated 
Curie temperate of the low temperature phase, and $0 
is a constant. 

The distortion d plays the same role as the relative 
change in volume (V — Vo)/Vq of the BR model (equation 
||TJ). In equation (JIJ Tq depends linearly on the change 
in volume, whereas now it depends quadratically on the 
distortion d, Tq — Tq(1 — Si d 2 ). This is the correct 
expansion of Tq, since the linear term in d disappears 
due to symmetry (+d and —d correspond to the same 
distortion). The terms in cr 2 , a 4 and cr 6 appear also in 
the BR model and lead to the first order transition at 
T p . The second order phase transition at T t = Td is 
generated by the d 2 and d A terms in the expansion. 

Also the variation of the magnetic moment with the 
distortion is contained in the model, although the authors 
find that the coefficient 82 is essentially zero. For Tq the 
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FIG. 5: Lattice parameters a and c (see figure [3J 

Curie temperature extrapolated from the high tempera- 
ture susceptibility above T t was considered (T = 285 K), 
and for Tb the temperature at which the inverse suscepti- 
bility has a minimum (To = 394 K). By construction this 
model yields the correct thermodynamic behavior, and 
also predicts the increase of the susceptibility between 
T p and T t by means of the reduction of the distortion 
with increasing temperature. 

Variations of this model give similar results. In 

reference [33J the full T-P-H phase diagram of MnAs is 
explained with a Landau-type expansion, where more or- 
der parameters are used. However also in this case a term 
equivalent to (T — To)rf 2 of equation is used in order 
to obtain the phase transition at Td- Despite the fact 
that Landau-type expansions give very good agreement 
between theory and experiments when the right parame- 
ters are used, they do not provide insights into the origin 
of the terms of the expansion, especially of the (T—To)d 2 
term. In reference |34| a basic justification of such term 
is given from first principles within a spin fluctuation 
theory constructed from a Hubbard Hamiltonian. It is 
shown that for a given volume the minimum of the free 
energy can lie at d = for a ferromagnetic state or at 
d for a paramagnetic state. 

In 1986-87 Motizuki and Katoh used spin fluctuation 
theory in order to explain the anomalous behavior of the 
susceptibility between T p and T t i 35 i 36 A Hubbard Hamil- 
tonian was used, with model density of states obtained 
from first principle calculations. They could qualitatively 
show that the susceptibility increases when going from 
the B31 structure to the B8i structure, again mainly due 
to the fact that Tc increases with increasing temperature. 

More recently various tight bindin g 37 ! 38 as well as first 
principlea 39 ! 40 ! 41 ! 42 ! 43 ! 44 ! 45 ! 46 ! 47 calculations have been 
performed for MnAs in the B8i structure. The results 
generally agree and compare well with the experiments. 
Only two studies on MnAs in the B31 structure are 
known to the authors i 4 ^ 4 ^ In reference 0| a descrip- 
tion of the paramagnetic state of the B31 structure is 




given by assuming that the paramagnetic state coincides 
with zero local magnetic moment of the Mn atoms. This 
does not correspond to the usual way of describing para- 
magnetism, which rather corresponds to constant local 
magnetic moments, whose directions in space change ran- 
domly in time due to spin fluctuations. 

The present work investigates the magnetic interac- 
tions across the various phase transitions of MnAs. An 
explanation of the magnetostructural properties in terms 
of first-principles calculations is given, thereby illustrat- 
ing the origin and providing a justification of the param- 
eters used by the different models. 



III. RESULTS 
A. DFT calculations 

First principles calculations within density functional 
theory (DFT) are performed using the pseudopotential 
code based on localized atomic orbitals SIESTA. 50 The 
generalized gradient approximation (GGA) is used for 
the exchange correlation potential, 5 i since it has been 
shown to give good structural properties for hexagonal 
MnAs. 39 i 41 i 42 In the valence we consider 4s, 4p and 3d 
orbitals for Mn, and 4s, 4p and Ad for As. For both 
Mn and As double zeta polarized local orbitals are used 
for the s and p angular momenta, whereas for the d or- 
bitals double zeta is used. The number of fc-points in the 
Brillouin zone is specified by a grid cutoff of 20 A. This 
corresponds to a 11x11x8 mesh for the B8i unit cell, 
giving approximately 1000 /c-points in the full Brillouin 
zone. For the B31 unit cell such cutoff yields a 8x11x7 
mesh. The real space mesh cutoff, which determines the 
density of the real space grid, is 300 Ry. 

After full relaxation of the unit cell to a pressure below 
0.1 kbar, and of the atomic positions to a force smaller 
than 0.01 eV/A the B31 unit cell in the ferromagnetic 
configuration relaxes to a B8i structure with a = 3.72 
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FIG. 6: B31 unit cell where only the four Mn atoms are shown 
displaced out of the B81 symmetry positions. 

A and c = 5.58 A. The experimental values at room tem- 
perature are a — 3.724 A and c = 5.707 A. Therefore 
the relative error is below 1% for a and —2% for c. The 
lattice parameters at 110 K can be extracted from figure 
[5] and arc a = 3.733 A and c = 5.677 A. This demon- 
strates that GGA-DFT reproduces rather well the zero 
temperature ground state. 

For a fixed c/a ratio of 1.54, which is close to the ex- 
perimental value at the first order phase transition, the 
energy is minimized for a = 3.695 A, which compares 
well with the results of other ab initio calculationsi^iiil^ 
Also the band structure and the density of states are sim- 
ilar to previous calculations. The magnetic moment per 
Mn atom is 3.4 /ib, and compares well with the measured 
value of 3.4 /xb^ 

The unit cell of the B31 structure contains 4 Mn 
atoms, allowing for 3 possible independent antiferro- 
magnetic configurations of the local moments of the 
Mn atoms. The different antiferromagnetic states are 

+ H , H h— and H h As a matter of notation 

+ H means that the atoms 1 and 2 in the unit cell 

have opposite magnetic moment than that of atoms 3 and 
4. The indices of the Mn atoms in the unit cell are de- 
fined in figure El A cell relaxation is performed for those 
three antiferromagnetic configurations. Table HT1 lists the 
obtained relaxed structures together with the total en- 
ergies per Mn atom as compared to the ferromagnetic 
ground state energy (E — -Efm)- The structure remains 
of the B8i type if the local moments are ferromagneti- 
cally aligned in the hexagonal plane, whereas it changes 
to the B31 type if the moments are antiferromagnetically 
aligned, with u (v) of the order of 0.05 b (c). There is 
also a slight displacement of the Mn atoms along c of at 
most 0.01 c, and of the As atoms along b of at most 0.01 
b. 

Generally it can be observed that the in-plane lattice 
parameters contract and the c-parameter expands for the 



antiferromagnetic states, resulting in a reduction of the 
volume V. The calculated values are similar to those 
given in reference [isj. 

The total energy for the ferromagnetic alignment is the 

lowest, although the H h configuration is higher of 

only 17 meV/Mn. This indicates that the system should 
evolve to one of the antiferromagnetic states under pres- 
sure, since those have a much smaller volume but only a 
slightly higher energy. 

Table|n]gi ve s also the distances between the first three 
nearest neighbor Mn atoms in the hexagonal plane e/2,1, 
d,2.i and ^2,3 (see figure 0}. Note that e?2,2 = a and it is 
not given explicitly. While these distances are all equal 
in the hexagonal case, they differ of as much as 1 A in the 
B31 structure. Large changes in the distance between the 
Mn atoms are possible since the nearest neighbor Mn- 
Mn separation in MnAs is well above the inter-atomic 
distance 2.61 A of bulk Mn^S which can be regarded as 
the minimal possible distance between Mn atoms. The 
distance between nearest neighbor Mn and As atoms lies 
between 2.46 A and 2.62 A for all the different configura- 
tions, and changes therefore much less than the Mn-Mn 
distance. 

The local magnetic moment on the Mn (/XMn) and 
As (^As) atoms, calculated using the atomic Mullikcn 
population, 53 is also given in table |n] The local moment 
on the Mn atoms ranges between 3.43 /ie for the ferro- 
magnetic configuration to 3.01 /zb for the H 1 — con- 
figuration. This reduction in the local moment is mainly 
due to the decrease of the cell volume, and the consequent 
increase of the hybridization between the Mn-d and As-p 
orbitals. 

In summary these calculations show that the distortion 
to the B31 structure is caused by an antiferromagnetic 
alignment of the local magnetic moments in the hexago- 
nal plane. 

B. Fit to Heisenberg energy 

In order to extract the various exchange parameters, 
calculations are performed for three different B31 super- 
cells in different local magnetic configurations. These 
supercells contain 8 Mn atoms, and are obtained by dou- 
bling the B31 unit cell along the a lattice vector (supercell 
1), along b (supercell 2) and along c (supercell 3). The 
calculated total energies are then fitted to a model Hamil- 
tonian. The most general form of Hamiltonian able to fit 
all possible energies for a system of N s collinear magnetic 
moments is 

E Su s 2 ,— — Eo~ J]] Tjy s ji s j2 ' ' ' s i" ^h-ji-,—,jv > 
"=1 ' h h iv 

(3) 

where Sj is the moment at site j (normalized to one), and 
the Js are coupling parameters. The Jjj,... parameters are 
symmetric under permutations of the indices, and each 
of the Jij.... is zero if any of the indices i and j are equal. 
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3.62 


6.29 


5.70 
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0.04 


3.12 


4.17 


3.33 


-0.03 
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TABLE II: Relaxed lattice parameters, local magnetic moment of the Mn and As atoms, and total energies per Mn atom for 
different spin configurations. 



The number of independent parameters is therefore 2 s . 
Since time reversal symmetry makes the system invariant 
under a global spin rotation only Js with an even number 
of indices are not zero and equation J2J reduces to 

^Si,32,... ^0 2 ^ ] Jij ^ ^ SiSj SfcSl Jijkl 

(4) 

I \^ T 

gj / J SiSj SkSlS m S n Jijklmn • • • 

i,j,k,l,m,n 

Here we neglect 4-moment coupling constants Jijki and 
higher thus reducing the model to an Ising-type. This ap- 
proximation neglects the dependence of the moment on 
each Mn atom as well as the small induced magnetic mo- 
ment of the As atoms on the orientation of the moments 
of the surrounding Mn atoms. Furthermore it is assumed 
that the coupling constants are independent from the an- 
gle between the magnetic moments, so that the fit can 
be extrapolated to vector magnetic moment to give a 
Heisenberg type energy 

E Sl , S2t ... = E - i s o s j J y- ( 5 ) 

In mean field theory the Curie temperature Tq for clas- 
sical Heisenberg exchanged magnetic moments is 

k B T c = "V 3 = J °A ( 6 ) 

j 

where fee is the Boltzmann constant. The quantity 
Jo = ■ Jo,j is the sum of the exchange coupling con- 
stants of a given magnetic moment with all the other mo- 
ments. In the following sections equation (jSJ) is always 
used to extract Curie temperatures, although it is well- 
known that the mean field approximation overestimates 
T C M. 

Calculations are performed for all the independent spin 
configurations of the supercell 1, and for a randomly cho- 
sen subset of those of supercells 2 and 3. 35 different 
configurations of the magnetic moments are used in to- 
tal. The energies are then fitted by a least-mean-square 
fit to the coupling parameters of equation JSJ. Since the 
system is metallic with the cZ-orbitals having finite den- 
sity of states at the Fermi level, the magnetic interaction 
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max 

FIG. 7: Variation of the various exchange quantities as func- 
tion of the number of coupling coefficients JV max included in 
the fit. (a) Standard deviation A of the energies resulting 
from equation JKJ as compared to the calculated DFT ener- 
gies per Mn atom, (b) Eo per Mn atom (equation ®), where 
the zero of energy is chosen as the energy of the ferromagnetic 
state, (c) Mean field Curie temperature Tc- 



is expected to have a long range character. For the cho- 
sen supercells it is possible to extract coupling constants 
up to the ninth nearest neighbor. The lattice parameters 
used are approximately those for ferromagnetic MnAs in 
the B8i structure at the phase transition temperature 
T p = 318 K (a=3.71 A, c/a=1.54). 

We have carefully tested the convergence of our results 
with the range of the Heisenberg exchange interaction. 
Figure [7| shows the standard deviation A of the energies 
resulting from equation JSJ as compared to the calculated 
DFT energies per Mn atom, the value of Eq per Mn atom 
and the mean field Curie temperature Tc as a function of 
the number of coupling coefficients iV max included in the 
fit. The standard deviation A decays monotonically, re- 
mains roughly constant for 7V max > 3, and then reaches a 
minimum value of around 5 mcV for -/V max = 9. This can 
be considered as the error resulting from neglecting high 
moment coupling constants. The value of Eq changes less 
over the whole range. 



8 



20 




FIG. 8: Exchange coupling parameters for JV max = 3 and (a) 

A^max = 9 as a function of the distance between the magnetic 

moments. 



Tq reaches a constant value of approximately 633 K 
for -/V max > 3. This indicates that the main contribu- 
tion arises from the first three nearest neighbor coupling 
constants. The experimental value of Tq for the low tem- 
perature phase ranges between T p = 318 K and T t = 400 
K. This means that our mean field Tq overestimates the 
experimental one by about a factor 2. 

Figure [S] shows the calculated exchange coupling con- 
stants as a function of the distance for two different fits 
counting respectively 3rd and 9th nearest neighbor cou- 
pling. The first three exchange constants remain nearly 
unchanged when going from 3rd to 9th nearest neighbor 
coupling. Interestingly the coupling parameters are pos- 
itive and therefore ferromagnetic up to d ~ 6. 5 A (fifth 
neighbor interaction). In what follows we consider only 
coupling parameters up to third nearest neighbors, as 
they give the main contribution to the properties of the 
material. 
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(b) 

FIG. 9: Schematic representation of the atomic positions of 
the Mn atoms together with the exchange constants for the 
B8i (a), and the B31 (b) structures. 



C. B8i to B31 distortion at T p 

The B8i to B31 structural phase transition at T p is 
investigated by calculating the Heisenberg coupling con- 
stants for different distorted cells. We start from the 
B8i with the experimental lattice parameters near T p 
(a = 3.71 A, b = \/~3 a, c — 1.54 a, u = v = 0) and 
then distort the cell linearly to the B31 structure. The 
amount of distortion d is given in percent, where d =0 
% stands for the lattice parameters of the ferromagnetic 
B8i cell at T p , and d = 100 % for the paramagnetic B31 
cell at T p (a = 3.676 A, b = 1.01%/3 a, c = 1.556 a, 
u = 0.02 b, v = 0.02 c). Calculations are done for dis- 
tortions between % and 220 %. Note that the volume 
decreases with increasing distortion. For these calcula- 
tions only the supcrcells 1 and 2 are used with a total of 
26 different spin configurations. The standard deviation 
of the fit is approximately constant for all the distortions 
and is of the order of 5 meV/Mn. 



Figure 9(a) shows the Mn atoms of the B8i structure 
coupled by first (Ji), second (Ja) and third (J3) nearest 
neighbor interaction. In the distorted B31 structure the 
three coupling constants Ji, J2 and J3 are split into six 
different constants due to symmetry loss. While there 
is still only one J\ coupling, the in-plane J2 splits into 
three different coupling constants 0/2,1, ^2,2 and J2,3, cor- 
responding to different distances between the Mn atoms 
in the hexagonal plane (see figure @J. Moreover also the 
third nearest neighbor coupling J3 splits into three dif- 
ferent constants, although two of them are between Mn 
atoms separated by approximately the same distance at 
T p , and so they are assumed to be identical. Hence J3 
effectively splits only into ^1 and J3,2- 

Figure^|shows the calculated values for the exchange 
parameters as a function of the distortion. For 0% dis- 
tortion the values of Ji,i, ^2,2 and ^3 are approximately 
equal reflecting the hexagonal symmetry. The values of 
Jz,i and 73,2 also should be identical although they differ 
by about 2 meV (note that in the fit we do not force the 
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B81 symmetry when determining the Js for the undis- 
torted structure). This can be assumed to be the error 
over the fit. Additional control fits were also performed 
for different subsets of the 26 spin configurations. The 
variation over the individual Js was of 20%, whereas the 
variation of Jo was always smaller than 6%. 

The value of J± remains approximately constant for 
all the distortions, reflecting the fact that the distance 
between the quasi-hexagonal layers remains roughly con- 
stant. In contrast the in-plane Js change and eventually 
become antiferromagnetic. In particular the coupling be- 
comes strongly antiferromagnetic for J2.1, i.e. for those 
Mn atoms that get closer in the hexagonal plane under 
distortion. Also J2,s, which couples the Mn atoms in- 
creasing their separation, is reduced and becomes anti- 
ferromagnetic for large distortions. Finally the coupling 
parameters J^i and J^ t 2 have only minor changes, with 
J31 becoming weakly antiferromagnetic for large distor- 
tions. 

The evolution of the coupling constants with the dis- 
tortion explains why for the H 1 — and H h spin 

configurations (table ^J, where the spins are antiferro- 
magnetically aligned in the hexagonal plane, the lowest 
energy is found for the B31 structure. In fact the re- 
laxed structure for both spin configurations is similar to 
a distortion of about 200 %. At this distortion the in- 
plane coupling constants J24 and ^2,3 become antiferro- 
magnetic, resulting in a reduction of the total energy as 
compared to the B81 structure. 

Figure 11(a) shows the relative change of the mean 



field Curie temperature Tc(d)/Tc(0) for the ferromag- 
netic state. Tq decreases monotonically with increasing 
the distortion. For 100 % distortion (B31 structure at 
T p ) r c (100%)/T c (0) = 0.67, demonstrating that when 
the phase transition from the B81 to the B31 structure 
occurs, the system in the B31 cell is already paramag- 
netic with very little magnetic order. The experimental 
Curie temperature T^, xp for the hexagonal cell at T = T p 
is not known, since the structure changes. 

Figure ll(b)| shows the total energy per B31 unit cell 
as a function of the distortion in the ferromagnetic (FM) 

and in the H h antiferromagnetic configurations (AF). 

This latter is the antiferromagnetic configuration giving 
the lowest total energy at its minimum among all the 
ones calculated along the considered distortion. The fig- 
ure also shows the value of Eq, the energy of the param- 
agnetic state (see equation (0). The zero in the energy 
scale is the energy of the ferromagnetic state for d = %. 

The ferromagnetic state has its energy minimum for 
d = %, and increases parabolically for increasing dis- 
tortion. This means that the B81 structure is the one 
with lowest energy in the ferromagnetic state. In con- 
trast the competing antiferromagnetic configuration has 
a minimum for about 180 % distortion, where the energy 
is lower than the ferromagnetic phase. The ground state 
of the system is therefore expected to be ferromagnetic 
for distortions up to about 100 %, and then to become 
antiferromagnetic for d >180 %. For distortions in be- 
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FIG. 10: Evolution of the exchange coupling constants when 
distorting the unit cell linearly from the B81 structure to 
the B31 structure, d — % represents the B81 structure 
at T p = 318 K, d — 100 % represents the B31 structure at T p . 
A positive (negative) value of J means ferromagnetic (anti- 
ferromagnetic) coupling. 



tween the ground state is expected to be a non-collinear, 
canted spin structure. 

Eq has a very flat minimum for distorted cells, reflect- 
ing the fact that the total energy increases for the ferro- 
magnetic state, but decreases on average for the antifer- 
romagnetic states. The minimum is found to be at about 
d ~ 100 %, which corresponds indeed to the lattice pa- 
rameters of the paramagnetic state above the phase tran- 
sition. This explains the structural change from the B81 
to the B31 structure when the MnAs magnetic state goes 
from ferromagnetic to paramagnetic. The phase transi- 
tion at T p can be seen as a jump from the minimum of 
the FM curve to the minimum of the Eq curve. 



D. B31 to B81 distortion at T t 

For temperature in between T p and T t the MnAs crys- 
tal structure continuously changes from B31 to B81. As 
mentioned in section III Al the phase transition temper- 
ature T t is usually identified as the temperature where 
the susceptibility and the specific heat have a maximum. 
This is at about 398 K. However, the distortion disap- 
pears at slightly higher temperatures as pointed out in 
reference [25|. Therefore, since the exact temperature 
for this second order structural phase transition is not 
known exactly, we introduce an operative definition and 
assume that the distortion disappears at a temperature 
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FIG. 12: Distortion do for the minimum of the paramagnetic 
ground state as a function of the lattice constant a. The dots 
are calculated values, the dashed line is a fit with equation 
0. 
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FIG. 11: (a) Relative change of the mean field Curie temper- 
ature (Tc(d) - r c (0))/r c (0) for the ferromagnetic state, (b) 
Total energy for one B31 unit cell for the ferromagnetic config- 
uration (FM), for the H h antiferromagnetic configuration 

(AF), together with Eq (equation JSJ), as a function of the 
distortion d. d = % and d = 100 % represent respectively 
the B8i and the B31 phase at T p «318 K. 



T s , at which the slope of the in-plane lattice constant 
as function of temperature a(T) changes discontinuously 
(see figure El- According to figure the lattice constant 
at T s is a(T s ) = h 3.699 A, and [da/dT) Ta « 0. At 
T t the same lattice parameter is a (T t ) «3.697 A, so that 
the difference in a is very small. 

The main reason for the second order phase transition 
at high temperatures is related to the lattice thermal 
expansion. The idea is that upon volume expansion, the 
ground state of the paramagnetic phase moves towards 
the hexagonal structure. We verify this hypothesis by 
calculating the minimum of Eq (_E™ m ) along a distortion 
of the cell transforming B8i to B31. In the calculation 



the volume of the cell is kept constant and we repeat 
the calculation for different volumes. This allows us to 
evaluate both E™ m and the corresponding distortion as a 
function of the volume. Since for T > T p MnAs is always 
paramagnetic, then the minimum of Eq corresponds to 
the stable distortion do at a given volume. In practice the 
change in volume can be described simply by the change 
in the planar lattice constant a, since both b/a and c/a 
do not deviate much from their value at T p . Thus we 
always consider b = v3 a and c = 1.556 a and the phase 
transition is investigated as a function of a only. 

The equilibrium distortion do as a function of a is 
presented in figure 1121 Indeed the distortion decreases 
with volume and it disappears for a between a = 3.695 
A and a = 3.71 A. Moreover we find 100% distortion 
for a ~ 3.66 A. These values agree rather well with the 
experimental ones, where the distortion disappears at 
about a(T s ) = 3.699 A, and 100% distortion is found 
at a(T p ) = 3.673 A. 

In order to interpret these results consider that the 
distortion is symmetric for ±d (Eo(d) — Eo(—d)), and 
therefore Eq can be expanded in even powers of the dis- 
tortion Eq((T) = ro + ri d 2 +r2 d A . Here are parameters 
to fit to the DFT calculations. In particular note that 
ro corresponds to the energy of the paramagnetic phase 
when the crystal is undistorted, i.e. it has the hexagonal 
structure. In this way the minimum of the Eo(d) curve is 

obtained for do = J—i^- For small distortions the pa- 
rameters r± and ri can be further expanded around a as 
T\ = ri,i(o — 5) and r 2 = r 2 ,o + r 2 ,i {a~d) + r 2 a (a-a) 2 . 
The calculated values for the leading terms are a = 3.706 
A, n,i = 62.4 10" 3 meV/A and r 2)0 = 1-29 10" 7 meV, 
(d is given in percent). The equilibrium distortion do up 
to first order in a is then 

Ma) = 1\ — — 6(a - a), 7= Ja—^-, (7) 

where Q(x) is the Heaviside function. With the values 
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FIG. 13: Distortion do for the minimum of the paramagnetic 
ground state as a function of temperature calculated with 
equation Q (with 7 = 1184 and a = 3.699) using a(T) taken 
from figure [K] 



FIG. 14: Energy expansion coefficient ro and i?o,min as 
function of the lattice constant. The dashed line shows 
eo+ei (a— ao) 2 , the solid line shows eo+£i (a — ao) 2 — £2 (a — a) 2 
(see equation ©)• 



of a, r\ i and r 2j o given above 7 = 947 is obtained. The 
resulting distortion is presented in figure IT2"1 

Interestingly if we use the equation Q to fit the ex- 
perimentally determined distortions at a(T s ) = 3.699A 
(d = 100%) and o(T p ) = 3.673A (d = 0), we obtain 
7 = 1184 and a = 3.699, both in good agreement with 
our calculated values. This suggests that the main ef- 
fects of the distortion to the B31 structure arise from the 
atomic displacement from the symmetry positions, and 
that small changes of the ratio of the lattice vectors, ne- 
glected in our calculations, play only a secondary role. 
Using the values 7 = 1184 and a = 3.699 of the two pa- 
rameters the evolution of the distortion as a function of 
temperature T can be obtained by inserting the data for 
a(T) from figurc[S]in equation (JJJ. The result is shown in 
figure^H It agrees well with the experimental result (see 
figure 3 in reference L.0t). The main difference however 
is that in reference [lfj the distortion do becomes zero at 
398 K, whereas in our results this happens only at 450 
K. This is due to our choice of T s = 450 K, which by 
definition sets the temperature where the distortion dis- 
appears. Close to the phase transition temperature fluc- 
tuations play an important role, so that for very small 
distortions close to the phase transition the description 
may not be valid. 

By using the computed values of ro, ri and r 2 the 
minimum of Eg (E™ ln =Eo(d = do)) is calculated as a 
function of the volume of the unit cell, and it is shown 
in figure [21 together with ro . Recalling that ro is the 
energy of the paramagnetic hexagonal phase, it also can 
be expanded as function of the lattice constant ro — eo + 
ei (a — ao) 2 , where ao is the equilibrium lattice constant 
of the hexagonal phase. This, combined with equation 
O, gives an expression for the energy minimum as a 
function of the lattice constant a 

E™ in (a) = e + e 1 (a- a ) 2 - e a (a - a) 2 0(5 - a), 



C2 



4r 2 . 



7546 



meV 



(8) 



where eo, e-y and ao are to be fitted from the calculations 
of ro (figure I14fl . The fitted values are eo = 264 meV, 
ei = 15935 meV/A 2 and ao = 3.65 A. From the equation 
© the energy minimum is easily found 



a 1 



e 2 a - a 
€1 - e 2 a 



(9) 



and by using the calculated parameters we estimate 
&min = 3.60 A. Since a m in < ao , we derive the important 
result that the distortion allows the volume to be further 
reduced as compared to the hexagonal phase. Further- 
more the curvature of the energy as a function of a is 



(a) 



da 2 



= 2ei - 2e 2 6 (5 - a) 



(10) 



which is also reduced by a factor 2e 2 when the structure 
is distorted. 

The effect of the thermal expansion on the lattice pa- 
rameter can now be modeled as a temperature dependent 
position of ao, ao = ao(T). The change of the lattice 
constant with temperature for the distorted phase can 
therefore be written as 



da u 



<9a n 



dT 



ei - e 2 dT 



(11) 



Since e 2 is smaller than ei the lattice expands with tem- 
perature faster for the distorted phase than for the undis- 
torted phase. If the calculated values for ei and e 2 are 
used the ratio ei/(ei — e 2 ) is found to be 1.90, which 
agrees well with the value of 2.25 extracted from figure 
[S] Near the phase transition, i.e. where d ~ 0, phononic 
effects due to the different curvatures of the energy and 
fluctuations should be considered. It is especially inter- 
esting that the change of the lattice constant with the 
temperature goes to zero near the phase transition tem- 
perature. 
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In the same way as Eq also Tq can be expanded as a 
function of the lattice constant and of the distortion 



OMn 



T c (a,d) = T c (a) 



l + K„ 



K d d 2 



(12) 



where K v and K d are parameters, and Tq (a) is the Curie 
temperature for the cell with lattice parameter a and 
where the atoms are in the hexagonal positions. The 
parameters are fitted by least mean squares to the cal- 
culated values of Tq obtained for six lattice constants 
ranging between 3.656 A and 3.78 A and for different 
distortions. We obtain Tc(a)=573 K, K v = 6.80 and 
K d = 2.62 10 -5 . At d = 100 % the relative change of the 
Curie temperature is T c (a, d = 100%)/T c (o) = 0.74 and 



corresponds roughly to the value of figure 11(a) In this 
case it is slightly larger due to the fact that the volume 
is kep t constant, whereas for the calculations of figure 
11 (a)| it shrinks with increasing the distortion. 

Next we calculate the dependence of the magnetic mo- 
ment on the distortion and on the unit cell volume. The 
dependence is again expanded to lowest order in a and d 



(13) 



We now have different ways of extracting the magnetic 
moment of the Mn atoms from our DFT calculations. 
One possibility is to take the total moment of the cell 
for the ferromagnetic spin configuration and divide it by 
the number of Mn atoms. In this way however the small 
induced moments of the As atoms are subtracted from 
the moment on the Mn. A second possibility is to take 
the average Mulliken spin population for the Mn atoms. 
The advantage of this method is that also antiferromag- 
netic configurations can be used to determine the average 
moment, and the induced moments of the As atoms are 
accounted for. The drawback however is that Mulliken 
populations are somewhat arbitrary as they depend on 
the basis set. 

By setting d to 3.699 A, the values obtained using 
the cell moment are /i(a) = 3.28 /xb, 0^=3.28 and 
a ft d=1.15 10~ 6 . Similarly from the average Mulliken 
population over all the magnetic configurations we ob- 
tain /z(o)=3.42 fiB, a M =3.48, and « 0. These results 
are rather similar to each other. Our analysis shows that 
the change of the magnetic moment is mainly due to the 
change in volume, whereas the distortion has basically 
no effect. 



E. Small distortions of the B8i structure 

In this section the dependence of Tq on the individ- 
ual lattice parameters and on the distances between the 
atoms is investigated for the B8i structure. Our ap- 
proach is to distort the cell orthorhombically but to leave 
the atoms in their high symmetry positions. Apart from 
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FIG. 15: Diagram of the positions of Mn and As atoms in 
orthorhombic MnAs. 



a general understanding of the phase diagram of MnAs 
this analysis is useful for predicting the behavior of MnAs 
when grown on a substrate. For instance when grown on 
GaAs(OOl) the substrate induces strain in MnAs, and 
the unit cell is slightly orthorhombically distorted^ This 
distortion does not correspond to the orthorhombic B31 
structure, since the atoms do not move out of the high 
symmetry positions. Moreover different growth orien- 
tations are possible, and the Curie temperature varies 
according ly* In addition reference presents experi- 
mental results showing that the phase transition temper- 
ature T p changes when strain is applied to the MnAs film. 
In that article we have compared our theoretical predic- 
tions for the dependence of T p on the lattice distortion 
with the experimental findings. In this section we refine 
and expand our previous analysis. 

When the cell is orthorhombically distorted the first 
three nearest neighbor coupling constants split into five 
different constants, corresponding to the ones of the B31 
cell (figure and with the only exception that now 
J2.1 = J2.3- The number of total energy calculations 
for the fit of the coupling parameters is 16. The change 
of Tq for each different distortion can be expressed as a 
function of the change of each single Mn-Mn and Mn-As 
distance in the unit cell. This gives 



5Jo 
Jo 



Sd u 



(14) 



where the sum goes over all five independent distances in 
the orthorhombic unit cell as defined in figure 1151 The 
dependence on the angles between the atoms is neglected. 

are the multiplicities of each distance <i M within one 
unit cell, and have the values iVMn.i = 4, Nmu,2 = 
4, A^Mn,3 = 8, N As ,i = 4 > ^As : 2 = 8. For the evalua- 
tion of the coefficients 21 different distortions are 
considered, including changes of volume, changes of the 
ratio of the different axes and different displacements of 
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the As atoms. For all the distortions the orthorhombic 
symmetry however was preserved. The best fit for Ka 
gives 



IV. DISCUSSION 



A. Curie temperature and susceptibility 



^Mn,l = 6.6, Kmtx,2 = 6.2, i^Mn,3 
^As,l = -7.5, K As .2 = -7.5. 



5.9, 



(15) 



The values of i^Mn.i an d Kmh.2 are almost identical, as 
one should expect from the symmetry and similarly for 
-Kas,i and K^ 2 - 

Equation (|14(l describes the fact that the change of 
Tq is the result of an interplay between the change of 
the Mn-Mn and Mn-As distances. The calculated 
show that, while an increase in the distance between Mn 
atoms increases Tq, an increase of the Mn-As distance 
decreases it. However note that the two distances can 
not be changed independently, thus the net change in Tq 
depends on the details of the distortion. 

Jo can also be expanded over the orthorhombic lattice 
parameters 



Jo 



3 

£ 

i=l 



Ki — -, with Ki 

a,i 



di dd„ 



where a\ = a, a 2 = b and 03 = c. The change of the 
distances between Mn and As atoms is not exactly known 
for the orthorhombic cell. However it is easy to show that 
to first order the position of the As atom in the cell does 
not influence Jq, since up to first order <Mas,2 = 2<Was,i 
when moving the As atom inside the cell. Therefore the 
As atoms can be assumed to remain in the high symmetry 
position. Assuming now K 2 — K3 and K4 = K$ (as 
imposed by symmetry), the general form of the Ki is 



K a = K h = 6K 2 + 
K c = AK X + 



96K 4 
16 + 3 c 2 ' 
36K 4 g 2 
16 + 3 5 2 : 



(17) 



where c = c/a. Using the average between K 2 and K3 
and c = 1.533 this gives 



K a = K h = A.9, K c = -l.O 



(18) 



Our results clearly show that stretching the unit cell 
along the basal plane raises T p (since K a > 0), while 
stretching along the c-axis lowers T p (K c < 0). An 
increase of the volume without distorting the cell results 
in an increase of the ferromagnetic exchange interactions 
and therefore of To, since K a is positive and larger in 
magnitude than K c . If the cell changes only its volume 
the expansion corresponds to the one of equation 112|) 
with a factor K v — 2K a + K c = 8.8. These results differ 
slightly from the ones given in reference [j| , since in that 
case we did not constrain K a to be equal to K b . 



It is now possible to analyze two peculiar properties of 
MnAs. The first is the anomalous behavior of the suscep- 
tibility x as a function of temperature between T p and 
Tf The second is the fact that although the Curie tem- 
perature for ferromagnetic MnAs has to be larger than 
T p , the Tc extrapolated from the susceptibility above T t 
is only 285 KiH In this section both these features are 
explained using the dependence of the susceptibility on 
the Curie temperature (x = x(Tc)), and the strong de- 
pendence of the Curie temperature on the lattice param- 
eters T c = T c {a,b,c,d). T C (T) and X (T) are therefore 
determined using the experimentally measured temper- 
ature dependence of the lattice vectors a(T),b(T) and 
c(T) and of the distortion d(T). This analysis also pro- 
vides a tool for extracting the parameters K v and K c i 
from experimental data. 

By generalizing the equations H12|) and (|16|) the Curie 
temperature Tc(a, b, c, d) can be written as 



(16) T G (a,b,c,d) 



l + K n 



a — a 6 — 6 



K r . C — C 



K d d 1 



(19) 



where we use the fact that K a = K b . We take the values 
for K a and K c from equation (fT£|l . K d = 2.62 10~ 5 and 
Tc,o = ?c(a) = 573 K are those calculated in section 
IIII DI The reference lattice parameters are chosen to be 
the lattice vectors at T s (d — 3.699 A, 6 = y/3 d,c = 
1.56 a). 

Similarly to the Curie temperature also the suscepti- 
bility is calculated in the mean field approximation. This 
is justified for T 3> Tc, a condition which is satisfied for 
paramagnetic MnAs. The molar susceptibility xm is then 
given by 



= 3 fc B s ( s + 1 )- 



(20) 



Na is the Avogadro's number, g as 2 is the Lande's factor 
for the free electron spin, is the Boltzmann constant 
and s is the atomic total spin. Note that the suscepti- 
bility has an additional implicit temperature dependence 
since Tc and s depend on the temperature through the 
lattice distortion. However in what follows we neglect 
the dependence of s on the lattice parameters so that Co 
is constant over all temperature. An analysis performed 
by relaxing this approximation gives similar results. 

Similarly to section lTlI Dl the model is further simplified 
by assuming that b/a and c/a are constant above T p . As 
indicated in equation JJJ d then is a function of the lattice 
constant, d(a) — 7-1/ (a — a) /a 0(a — a). Moreover we 
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have shown that the experimental distortion as function 
of the lattice constant is well reproduced when 7 = 1184. 
Therefore it is now possible to express Tc and Xm 1 as a 
function of the lattice constant a only 



Xm 1 ^) = 
1 



T - Tc.o 1 + {K v + K d 7 2 6 (5 - a)) 



(21) 



where K v = 2K a + K c = 8.8. This equation shows that 
if the lattice expands strongly with temperature Xm 1 ^e- 
creases. 

For temperatures between T p and about 390 K, as well 
as above T s , a increases approximately linearly with tem- 
perature (see figure [SJ) and can therefore be written as 



a(T) = o(T ) 



1 



T-Tn 



T 



(22) 



where the experimental values for the coefficients are 
T ,+ = T s = 452 K, a(T s ) = a = 3.699 A and a+ = 
0.0126 for temperatures above T s (the index "+" denotes 
the high temperature region above T s ), and To t - —T p — 
318 K, a(T p ) = 3.673 and a_ = 0.0284 for temperatures 
between T p and about 390 K (the index "— " denotes the 
intermediate temperature region) . By inserting equation 
(|22|l into equation l|21|) we obtain for the high tempera- 
ture region above T s 



Xm(T) = -J—{T - T c , e g), 



(23) 



with 



C e ff — 



Tc,ett — 



l-K v a+^A 
1 + K v a+ 



Co, 



(24) 



l-K va+ T -§A Tc '°- 



Tc iC ff and C e g are the experimentally accessible quanti- 
ties for the high temperature susceptibility, and due to 
the expansion of the lattice they are different from Tc,o 
and Co- The experimentally measured values are Tc je ff= 
285 K and C cff = 3.12 10~ 5 m 3 K, 14 which corresponds 
to an effective magnetic moment of 3.57 /ib- From C e ff 
and To, off the values of Tq.q = lb (2) and Co can now be 
obtained 



C n =- 



1 — K v a.\ 



1 - K v a+ 1 



T 3 



a 



eff ■ 



(25) 



Tc,o 



I - K„ou. 1 



C.eff- 



All the variables on the right hand side of equation l|25|) 
can be obtained from experiment except K v . For small 
K v the difference between Tc.o and Tc c ff is proportional 



to K v . Since T s is larger than the experimental value of 
7c. c ff the effect of the thermal expansion of the hexag- 
onal structure is a reduction of the slope of the inverse 
susceptibility as a function of temperature, as well as a 
reduction of the extrapolated Curie temperature as com- 
pared to the real Curie temperature. 

In the region where linear expansion holds the slope of 
the inverse susceptibility above T t is 



d X M ]+ 
dT 



Tc,o 



(26) 



whereas for in the intermediate temperature region above 
is 

^Xm 1 - 1 f, , T , ?\ a o T c ,o 



T p it is 



dT 



(27) 



In both regions there is a reduction of the slope due to 
the expansion of the lattice. However the reduction is 
much larger for xm,- than for xm,+ , since there is the 
additional term proportional to Kd due to the distortion, 
and also a_ > a + . As a rough approximation it can 

d x ~ 1 _ 

be assumed that (ao Tc,o)/(a To,-) ~ 1, so that 
becomes negative for 

(K v + K d 7 2 ) a_>l. 



dT 



(28) 



The values of a and 7 are determined experimentally 
and describe how the structure changes with tempera- 
ture, whereas K v and Kd describe how Tc varies for dis- 
torted cells. By using our calculated values for K v and 
Kd we obtain (K v + Kd 7 2 ) a- = 1.29. This is indeed 
larger than one. Therefore we do predict a negative slope 
for the inverse susceptibility in the intermediate temper- 
ature region. 

Finally we extract the values for K v and Kd from the 
experimental behavior of the Curie temperature. Since 
the ratio between K a and K c can not be obtained from 
the thermal properties of MnAs, it is therefore assumed 
that K c /K a = —1/4.9 = —0.2 is fixed and corresponds 
to our calculated value. For the hexagonal cell (d = 0) 
the equation (|19fl reads 



T c (a,c) = T, 



c,o 



l + K, k, 



(29) 



where k a = 2/(2 + K c /K a ) = 1.11 and k c = 
(K c /K a )/(2 + K c /K a ) = -0.11. This has to be valid 
for all temperatures where the cell is hexagonal, there- 
fore it can not be assumed that the ratio between c and a 
is constant since it changes abruptly from 1.533 to 1.556 
at T p . By using the expression for Tc.o from equation 
(|25H we obtain K v 



K.„ 



Tc(a,c) 

Ifa, off 



T c (a,c) 



Tc(a,c) 



+ k a 



(30) 
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As reference Curie temperature the extrapolated value 
to room temperature is used, which can be estimated 
to be about T c (a = 3.724 A, c = 1.533a)=360 K. By 
inserting the experimental values for the parameters on 
the right hand side of the equation 130|) we obtain K v — 
18.1. This is about twice as big as our predicted value. 
The disagreement may be due to the fact that the ratio 
between K c and K a has been fixed for our calculated 
value. By inserting this value for K v in equation i|25|) 
we obtain Tq $ = 311 K, the Curie temperature for the 
lattice parameters at T = T s . 

In order to extract Kd from experiments we use the 
relative change in the slope of the inverse susceptibility 
around Tt. This is, according to the equations l|26|l and 



dT 



dXM 1 - 
dT 



1 - K v 



an Th.o 



dT 



1 



(31) 

and increases with increasing K v and K4. All the vari- 
ables in this equation can be derived from the experi- 
mental measurements, except K v and Kd- Experimen- 
tally different values are found for the relative change 
of the slope (left hand side of equation l{3*Tj) )i 25 i 55 i 56 i 57 
These are all of the order of 1.44. By using this value for 
the relative change of the slope and the previously calcu- 
lated value K v = 18.1, K d is found to be 1.78 10" 5 . This 
value agrees approximately with our predicted value of 
2.62 1(T 5 . 

In conclusion the table II I II summarizes the parame- 
ters calculated in this work, by comparing our ab ini- 
tio results obtained from the DFT calculations and the 
Heisenberg model, with the results obtained by fitting to 
the experimental data. In general there is a fair agree- 
ment, although the the DFT results underestimate K v 
and overestimate Kd- 





Tc,o 


Ki 


K v 


Ka 




7 


a (A) 




(K) 


(io- 


5 ) 










DFT 


579 


2.62 


8.8 


4.9 


-1.0 


947 


3.706 


FIT 


311 


1.78 


18.1 


10.1 


-2.1 


1184 


3.699 



TABLE III: Main parameters used in the description of the 
phase diagram of MnAs. We compared results obtained from 
ab initio calculations and Heisenberg model (DFT), with 
those of the best fit to the experimental properties (FIT) . 



B. Phase diagram 

It is now possible to draw a qualitative description of 
the phase diagram of MnAs (see figure EJ) ■ Since the 
magnetic interactions are strongly dependent on the vol- 
ume, the effects of the thermal expansion of the lattice 



are important. First the state of MnAs at zero pres- 
sure is considered. For P = and temperatures below 
155 K both ferromagnetic and canted spin structures can 
be stabilized. The ground state is the ferromagnetic B8i 
structure, however also the B31 structure is stable, be- 
cause at these low temperatures the thermal energy is not 
large enough to induce a phase transition from a canted 
spin structure to a ferromagnetic one. As the tempera- 
ture increases above 155 K the thermal fluctuation may 
induce the phase transition, and therefore only the ferro- 
magnetic state survives. 

As the temperature is further increased and the system 
becomes paramagnetic, the B31 structure becomes the 
stable state, as described in section ITlI CI (figure 11(b) I. 



MnAs then remains paramagnetic for all temperatures, 
but because of the thermal expansion the orthogonal dis- 
tortion decreases, until it vanishes when a critical volume 
is reached (see figure lT2*|) . Then the volume becomes large 
enough to stabilize the B8i structure even for a para- 
magnetic state, leading to the second structural phase 
transition at T t . 

As pressure is applied to the system the volume is 
reduced and thus the ferromagnetic exchange interac- 
tion decreases. Therefore T p decreases with pressure, 
while Tt increases. At a temperature of approximately 
230 K the ferromagnetic state is stable up to a pres- 
sure of 2kbar. For high pressures above 2.5-4 kbar the 
volume is small enough that canted spin structures be- 
come the ground state. The temperature at which these 
states become paramagnetic is approximately 230 K. It 
increases slightly with pressure, reflecting the fact that 
the antiferromagnetic interaction increases as the volume 
is reduced. 



V. CONCLUSIONS 

We have investigated by means of ab initio electronic 
structure calculations the magneto-structural properties 
of MnAs. The stable structure for the ferromagnetic state 
is found to be the B8i structure. However if antiferro- 
magnetic alignment in the hexagonal plane is imposed 
the B31 structure becomes more stable. By fitting the 
DFT total energies of different magnetic configurations to 
a Heisenberg type energy it is shown that the main con- 
tributions to the physical properties originate from the 
exchange coupling parameters up to 3rd nearest neigh- 
bor. The Curie temperature was calculated in the mean 
field approximation, with values approximately twice as 
large as the experimental ones. 

The main assumption of the phenomenological model 
of Bean and Rodbellii that the ferromagnetic exchange 
coupling parameters increase when the volume is in- 
creased has been confirmed (equation Q18|0 using this 
analysis. However it has been shown that the exchange 
interactions depend not only on the volume, but that 
the orthogonal distortion to the B31 structure plays an 
important role. For the experimentally observed distor- 
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tions some of the in-plane exchange coupling coefficients 
become antifcrromagnctic. This is the reason for the sta- 
bility of the B31 structure for those configurations of the 
magnetic moments that have an antiferromagnetic com- 
ponent in the hexagonal plane. Different canted spin 
structures are expected to minimize the energy for dif- 
ferent distortions, since there are both positive and neg- 
ative exchange coefficients depending on the amount of 
distortion of the B31 structure. 

Furthermore it has been shown that for paramagnetic 
states the B31 structure is stable at small volumes, while 
the B81 structure is stable above a critical lattice con- 
stant of about 3.7 A. This explains the second order phase 
transition at T t , since at that temperature the lattice con- 
stant crosses this critical value. 

The Curie temperature has been expanded as a func- 
tion of the lattice vectors and of the amount of distor- 
tion. An increase in the volume leads to an enhancement 
of the Curie temperature, while an increase of the dis- 
tortion leads to a reduction. With these results in hand 
the increase of the susceptibility between T p and T t has 



been explained as the result of the increase of the Curie 
temperature due to the change of the structure from the 
B31 to the B81 and to the increase of the volume. By us- 
ing the experimental variation of the lattice parameters 
with increasing temperature the susceptibility is indeed 
found to increase between T p and T t . 

A fit of the dependence of the Curie temperature on 
the lattice parameters to best reproduce the experimen- 
tal behavior is also given. The calculated values agree 
within a factor of two with the values obtained by the 
ab initio calculations. Our results are in agreement with 
the phenomenological models based on the Bean Rodbcll 
idea. But now the various parameters used are derived 
from first principles and therefore validated. 
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